Noise activated switching in a driven, nonlinear micromechanical oscillator 
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We study noise induced switching in systems far from equilibrium by using an underdamped 
micromechanical torsional oscillator driven into the nonlinear regime. Within a certain range of 
driving frequencies, the oscillator possesses two stable dynamical states with different oscillation 
amplitudes. We induce the oscillator to escape from one dynamical state into the other by intro- 
ducing noise in the excitation. By measuring the rate of random transitions as a function of noise 
intensity, we deduce the activation energy as a function of frequency detuning. Close to the criti- 
cal point, the activation energy is expected to display system-independent scaling. The measured 
critical exponent is in good agreement with variational calculations and asymptotic scaling theory. 

PACS numbers: 05.40.-a, 05.40.Ca, 05.45.-a, 89. 75. Da 
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Fluctuation-induced escape from a metastable state is 
an important problem that is relevant to many phenom- 
ena, such as protein folding and nucleation in phase tran- 
sitions. For systems in thermal equilibrium, the escape 
rate can be deduced from the height of the free-energy 
barrier Q . The barrier decreases as the control parame- 
ter fx approaches a critical (bifurcational) value \i c where 
the metastable state disappears. It has been established 
theoretically and experimentally 0, El that, in the sim- 
plest and arguably most common case of the saddle-node 
(spinodal) bifurcation 0, the barrier height scales as 
(fx — ^ c ) 3 / 2 - Much less is known about escape in sys- 
tems far from thermal equilibrium 0, 0, ■ Such sys- 
tems are not characterized by free energy, and the scaling 
behavior of the escape rate near a saddle-node bifurca- 
tion has not been studied experimentally. Recently the 
problem of escape far from equilibrium has attracted sig- 
nificant experimental attention in the context of systems 
where multistability itself arises as a result of strong pe- 
riodic modulation. Escape was studied in parametrically 
driven electrons in a Penning trap [jj, doubly clamped 
nanomechanical oscillators Qfl lll| and radio frequency 
driven Josephson junctions |12j . 

We report here our investigation of noise-activated 
switching in systems far from equilibrium. By using 
a well-characterized system, an underdamped microme- 
chanical torsional oscillator periodically driven into non- 
linear oscillations, we study the dependence of the escape 
rate on the control parameter as it approaches the criti- 
cal value and reveal the scaling of the activation energy 
of escape in a system far from thermal equilibrium. The 
strongly-driven micromechanical oscillator has two sta- 
ble dynamical states with different oscillation amplitude 
within a certain range of driving frequencies. We induce 
the oscillator to escape from one state into the other by 
injecting noise in the driving force. By measuring the 
rate of random transitions as a function of noise inten- 
sity, we demonstrate the activated behavior for switch- 
ing and deduce the activation energy as a function of 
frequency detuning. Close to the bifurcation frequency 
where the high-amplitude state disappears, the activa- 
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FIG. 1: (a) A cross-sectional schematic of the micromechan- 
ical torsional oscillator with electrical connections and mea- 
surement circuitry (not to scale), (b) Normalized frequency 
response of sample A for excitation voltages amplitudes of 
9.1 ijV (squares) and 86.5 /J.V (circles). The dotted line repre- 
sents a fit to the data at smaller excitation using the response 
of a damped harmonic oscillator. For the large excitation, two 
dynamical states coexist from 3278.76 Hz to 3280.8 Hz. The 
dashed line fits the data to a damped oscillator with cubic 
non-linearity [3, yielding (3 = -1.2 x 10 14 rad~ 2 s~ 2 . 



tion energy is predicted by variational calculations and 
asymptotic scaling theory to display system-independent 
scaling 0, ITol Hq . Our measurement of activation en- 
ergy as a function of frequency detuning yields critical 
exponents of 1.38 ± 0.15 and 1.40 ± 0.15 respectively for 
two samples with different parameters, in good agree- 
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ment with the theoretical predictions for escape in non- 
equilibrium systems near a spinodal point. 

In our experiment, measurements were performed on 
two micromechanical oscillators (samples A and B). 
Each oscillator consists of a moveable polysilicon plate 
(500 fiui by 500 /im by 3.5 fim) that is supported by 
two torsional rods (40 \xm by 4 /im by 2 /im). The 
other ends of the torsional rods are anchored to the sili- 
con substrate. A 2-/im-thick sacrificial silicon oxide layer 
beneath the plate was etched away, leaving the plate 
suspended between the torsional rods. The moment of 
inertia I for the plate is 4.3 x 10~ 17 kg to 2 , while the 
torsional spring constants are 1.83 x 10~ 6 N m rad -1 
and 1.85 x 10~ 6 TV to rad -1 for samples A and B re- 
spectively. Beneath the top plate, there are two fixed 
electrodes (500 [im by 250 fim) on each side of the tor- 
sional springs. One of the electrodes is used for exciting 
the torsional oscillations while the other electrode is used 
for detecting the oscillations. More details about the os- 
cillators can be found in Refs. 0,0]. 

Figure QJi shows a cross-sectional schematic of the os- 
cillator with electrical connections and measurement cir- 
cuitry. The application of a periodic voltage with DC 
bias Vdd to one of the electrodes leads to an electrostatic 
attraction between the grounded top plate and the elec- 
trode. Torsional oscillations of the top plate are excited 
by the periodic component of the electrostatic torque. 
The detection electrode is connected to a dc voltage Vd C 2 
through a resistor R. As the plate oscillates, the ca- 
pacitance between the plate and the detection electrode 
changes. The detection electrode is connected to a charge 
sensitive preamplifier followed by a lock-in amplifier that 
measures the signal at the excitation frequency. Mea- 
surements for both samples were performed at pressure 
of less than 2 x 10~ 7 torr. Sample A was measured at 
liquid nitrogen temperature and sample B was measured 
at liquid helium temperature. The quality factors Q are 
about 2, 000 and 8, 000 for samples A and B respectively 
at the corresponding measurement temperatures. 

The strong nonlinearity in our torsional oscillator de- 
rives mainly from the distance dependent electrostatic 
interaction between the top plate and the excitation elec- 
trode. The equation of motion of the oscillator is given 
by 0: 

6 + 2^ + ^6 = j, (1) 

where 6 is the angular rotation of the top plate, 7 is the 
damping coefficient, ujq is the natural frequency of the 
oscillator and r is the driving torque. The driving torque 
arises from the electrostatic interaction between the top 
plate and the driving electrode. For small rotation angle 
9, t can be written as: 

r = bF, (2) 

where b (~ 125 ^to) is half the distance from the axis of 
rotation to the edge of the top plate. The electrostatic 



force F is evaluated with separation d — b0: 

p — ^ AV2 (3) 

2{d-bey () 

As shown in Fig. QJl, d (~ 2 /dm) is the separation be- 
tween the top plate and the electrode in the absence of 
external torque, A is the area of the electrode and eo is 
the permittivity of free space. The excitation voltage V 
consists of three components: 

V = V dcl + V ac sin{ujt) + V noise (t). (4) 

The three terms on the right side of Eq. (0J represent the 
dc voltage, periodic ac voltage with angular frequency 
lo and random noise voltage respectively. Vdcl is chosen 
to be much larger than V ac and V no i se to linearize the 
dependence of F on V ac and V no i se . The strong distance 
dependence of the electrostatic force leads to nonlinear 
contributions to the restoring torque. A Taylor expansion 
of the electrostatic force F about d in Eqs. JQ) and (J2J) 
leads to: 

6+2 1 6+[Lol-ri]6+a6 2 +l36 z -C = Esin(ut)+N(t), (5) 

wWp™— 3b 3 c A v 2 a_ 2b 4 epA v 2 _ b 2 e A v 2 

C = Ii¥ V dci' E = be °2J^ r1 V ^ is the effective amplitude 
of the periodic excitation and N(t) = be °^i cl V noise (t) is 
the effective noise in the excitation. The linear term in 
the Taylor expansion of F produces a shift in the nat- 
ural frequency of oscillation, while the nonlinear effects 
arising from a and [3 are characterized by a constant 
[l3j | k = 3/3/8wo — 5a 2 /12o;o 2 - For our device geometry, 
~/3/a ss 330 and /3/uj 2 w 10000. As a result, the contri- 
bution of the quadratic nonlinearity (the a term) to k 
is negligible compared to the cubic nonlinearity (the (3 
term) and our device can be regarded as a Duffing oscil- 
lator in the absence of injected noise. 

First, we focus on the response of the oscillator with 
no injected noise in the excitation. We show in Fig. ^ 
the frequency response of sample A at two different ex- 
citation amplitudes. Both responses have been normal- 
ized so their respective maximum amplitudes are one. 
The squares represent the response of the oscillator at 
the smaller excitation. The resonance peak is fitted well 
by the dotted line that corresponds to the response of a 
damped harmonic oscillator. As the periodic excitation is 
increased, the cubic term in Eq. JSJ leads to nonlinear be- 
havior in the oscillations. The resonance curve becomes 
asymmetric with the peak shifting to lower frequencies, 
consistent with a negative value of (5. At a high enough 
excitation, hysteresis occurs in the frequency response, as 
shown by the circles in Fig.^. Within a certain range of 
driving frequencies, there are two stable dynamic states 
with different oscillation amplitude and phase. Depend- 
ing on the history of the oscillator, the system resides in 
either the high-amplitude state or low-amplitude state. 
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In the absence of fluctuations, the oscillator remains in 
one of the stable states indefinitely. 

When sufficient noise is applied in the excitation, the 
oscillator is induced to escape from one dynamic state 
into the other. Since this driven, bistable system is far 
from thermal equilibrium and cannot be characterized by 
free energy, calculation of the escape rate is a non-trivial 
problem. Theoretical analysis 0, 0] suggests that the 
rate of escape T at a particular driving frequency depends 
exponentially on the ratio of an activation energy E a to 
the noise intensity In- 



fcce^. 
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Close to the bifurcation frequency, where the high- 
amplitude state disappears, the activation energy is ex- 
pected to display system-independent scaling: 



E a oc (y - v h y oc (Ai/)', 



(7) 



where the detuning frequency Av is the difference be- 
tween the driving frequency v and the bifurcation fre- 
quency ff,. The activation energy is predicted GHH to 
increase with frequency detuning with critical exponent 
p = 3/2. To our knowledge, this scaling relation has not 
yet been observed, even though it is expected to occur 
in a number of non-equilibrium systems. We describe 
below our comprehensive experimental investigation of 
activated switching from the high-amplitude to the low- 
amplitude state for two micromechanical oscillators with 
different resonant frequencies and damping coefficients. 
The critical exponents measured for both samples were 
in good agreement with theory. 

In our experiment, we induce transitions from the high- 
amplitude state to the low-amplitude state by inject- 
ing noise in the excitation with a bandwidth of 100 Hz 
centered about the resonant frequency. The bandwidth 
of the noise is much larger than the width of the res- 
onance peak. We chose a sufficiently large sinusoidal 
excitation so that the hysteresis loop exceeds twice the 
resonance peak width. Figure |2t shows typical switch- 
ing events at an excitation frequency of 3278.81 Hz 
for sample A where the oscillator resides in the high- 
amplitude state for various durations before escaping to 
the low-amplitude state. Due to the random nature of 
the transitions, a large number of switching events must 
be recorded to determine the transition rate accurately. 
During the time interval between switching events in Fig. 
Et; the oscillator is reset to the high-amplitude state us- 
ing the following procedure. First, the noise is turned off 
and the driving frequency is increased beyond the range 
of frequencies where bistability occurs (> 3280.8 Hz as 
shown in Fig. \l]p). The driving frequency is then de- 
creased slowly towards the target frequency so that the 
oscillator remains in the high-amplitude state. Once the 
target frequency is reached, the noise is turned back on 
and the time for the oscillator to escape from the high- 



amplitude state is recorded. This process is then re- 
peated multiple times to accumulate the statistics for 
switching. Such a procedure is necessary because the en- 
ergy barrier for transitions from the low-amplitude state 
back to the high-amplitude state is much larger than the 
barrier for transition in the opposite direction. Thus, 
noise induced transitions from the low-amplitude state 
to the high-amplitude state will fail to occur in the du- 
ration of the experiment and the oscillator must be re- 
set to the high-amplitude state using the steps described 
above. Figure^ shows a histogram of the residence time 
in the high-amplitude state before a transition occurs. 
The exponential dependence on the residence time indi- 
cates that the transitions are random and follow Poisson 
statistics as expected. 
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FIG. 2: (a) In the presence of noise in the excitation, the 
oscillator switches from the high-amplitude state to the low- 
amplitude state at different time intervals. The system is 
reset to the high-amplitude state between switching events. 
The detuning frequency is Av — 0.05 Hz. (b) Histogram of 
the residence time in the upper state before switching occurs, 
at a detuning frequency of 0.05 Hz for sample A. The dotted 
line is an exponential fit. 



To determine the activation energy for a particular de- 
tuning frequency, we record a large number of transitions 
for multiple noise intensities (In)- The average residence 
time at each noise intensity is extracted from the ex- 
ponential fit to the corresponding histograms. Figure |3| 
plots the logarithm of the average transition rate as a 
function of inverse noise intensity. The transition rate 
varies exponentially with inverse noise intensity, demon- 
strating that escape from the high-amplitude state is ac- 
tivated in nature. According to Eq. ©, the slope in Fig. 
|3lyields the activation energy for escaping from the high- 
amplitude state at the particular detuning frequency. 

We repeat the above procedure to determine the ac- 
tivation energy for other detuning frequencies Av (Av 
is the difference between the driving frequency and the 
bifurcation frequency at which the high-amplitude state 
disappears) and show the result in Fig. 0Ji for sample A 
and Fig. 0Jd for sample B. All the detuning frequencies 
chosen for sample A are smaller than its resonance peak 
width while the maximum detuning frequency for sample 
B is about 4 times its resonance peak width. Fitting the 
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FIG. 3: Logarithm of the transition rate from the high- 
amplitude state as a function of inverse noise intensity at a 
detuning frequency Av of 0.25 Hz for sample A. The slope of 
the linear fit yields the activation energy. 



activation energies with a power law dependence on the 
detuning frequency yields critical exponents of 1.38±0.15 
for sample A and 1.40 ± 0.15 for sample B respectively. 
Despite the different resonant frequencies and a factor of 
4 difference in damping, the critical exponents obtained 
for both samples are in good agreement with theoretical 
predictions [lH ITsL Hij. To our knowledge, this repre- 
sents the first experimental verification of universal scal- 
ing for activated escape close to saddle-node bifurcations 
in systems far from thermal equilibrium. 
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FIG. 4: Dependence of the activation energy on detuning 
frequency for (a) sample A and (b) sample B. The solid lines 
are power law fits, yielding critical exponents of 1.38 ± 0.15 
and 1.4 ± 0.15 respectively. 

We remark that recently Aldridge and Cleland 01 
measured noise-induced switching between dynamical 
states in a nanomechanical beam. They found a 
quadratic dependence of the activation barrier on the 
distance to the critical point where the two stable states 
of forced vibrations and the unstable periodic state all 
merge together. Such quadratic dependence arises [l5]| 
when the parameters are changed along the line where 
the populations of the two stable states are equal to each 
other. This requires changing simultaneously both the 



amplitude and the frequency of the driving field. In con- 
trast, in our experiment, we approach a bifurcation point 
where a stable large-amplitude state and the unstable 
state merge together, while the stable small-amplitude 
state is far away. We vary only one parameter, the de- 
tuning frequency, while maintaining the periodic driving 
amplitude constant. We found that the activation bar- 
rier for escape is reduced to zero with critical exponent 
of 3/2. 

In conclusion, we demonstrated the activated behavior 
of noise induced switching of a thermally nonequilibrium 
system, a nonlinear underdamped micromechanical tor- 
sional oscillators modulated by a strong resonant field. 
The measured critical exponent for the activation energy 
near the bifurcation point agrees well with the predicted 
value of 3/2, verifying the system-independent scaling of 
the activation energy in the vicinity of the bifurcation 
point. Other systems that are far from equilibrium, in- 
cluding RF-driven josephson junctions |l2lll9| . nanomag- 
nets driven by polarized curre nt ] 20l| and double barrier 
resonant tunneling structures [21|, are expected to obey 
the same scaling relationship near the spinodal point. 

We thank M. I. Dykman and D. Ryvkinc for useful 
discussions. 
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